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Abstract. Predictions from Hartree-Fock (HF), Popov (P), Yukalov-Yukalova (YY) 
and t-matrix approximations regarding the thermodynamics from the normal to the 
BEC phase in weakly interacting Bose gases are considered. By analyzing the 
dependence of the chemical potential /i on temperature T and particle density p 
we show that none of them predicts a second-order phase transition as required by 
symmetry-breaking general considerations. In this work we find that the isothermal 
compressibility kt predicted by these theories does not diverge at criticality as expected 
in a true second-order phase transition. Moreover the isotherms ^ = fi{p,T) typically 
exhibit a non-singled valued behavior in the vicinity of the BEC transition, a feature 
forbidden by general thermodynamic principles. This behavior can be avoided if a 
first order phase transition is appealed. The facts described above show that although 
these mean field approximations give correct results near zero temperature they are 
endowed with thermodynamic anomalies in the vicinity of the BEC transition. We 
address the implications of these results in the interpretation of current experiments 
with ultracold trapped alkali gases. 
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1. Introduction 

In the recent past it has been widely discussed that Bose-Einstein Condensation (BEC) 
in interacting Bose systems, i.e, the passage from the normal to the condensed phase 
where all particles occupy a single-particle state, shows an spontaneous U{1) gauge 
symmetry breaking with the condensate fraction Nq/N playing the role of the order 
parameter [H [21 [Sj H]. The overwhelming task faced when trying to solve the full inter- 
acting quantum gas has motivated the search of physical approximations for the dilute 
and weakly interacting gas O |6] . This search has been also encouraged by the experi- 
mental realization of BEC in alkali gases[71[H] since BEC occurs in the regime of s wave 
scattering where the interaction potential U{fi — 7^2) may be approximated by a contact 
potential U{ri—r2) = Uo6^{ri—r2)^. The introduction of this simplification has lead to 
a formulation valid at low densities and temperatures near T = 0. This approximation 
formally known as the theory of weakly interacting Bose gases has become the standard 
tool for analyzing the BEC transitions for gases confined either in a box of volume V 
or trapped through external potentials in optical setups. 

In spite of the great success achieved, the theory of weakly interacting Bose gas still 
awaits for a formal and complete analytical solution. In order to extract some partial re- 
sults different approximations have been developed in an upper 'layer' build up on top of 
the theory of weakly interacting gases. This 'layer' contains additional approximations 
to the already approximated theory of weakly interacting gases [21 [ini [HI [121 [iHl [lH [I5] . 
In this work we focus our attention to four of the approximations most widely used to 
describe the interacting Bose gas in the dilute approximation. These are Hartree-Fock 
(HP), Popov (P), Yukalov-Yukalova (YY) and t-matrix approximations. As we show in 
this work none of these approximations exhibit BEC as a second-order phase transition. 
The fact that HP neither exhibits a second-order phase transition nor displays the cor- 
rect energy spectrum has been known for a while flU\ ITT] . Here, we show that Popov 
(P), t-matrix approach|21 [13] and the recently introduced Yukalov-Yukalova (YY) [12] 
approximations also fail to predict a second-order phase transition. In addition we show 
that all of these theories display an unstable region implying that the BEC transition 
is at best, a first-order phase transition. 

To achieve this goal we proceed as follows. We start by calculating the equation of 
state /i = /i(T, p), that is, the dependence of the chemical potential /i on temperature 
T and particle density p = N/V, with N the number of particles and V the volume 
of the system in each of the approximations considered. Equilibrium thermodynamics 
ensures that the isothermal compressibility kj- can be obtained from the equation of 
state /i = n{T, p) through the relationship = p'^{dfi/dp)T from which information 
regarding the order of the transition predicted by these theories can be inferred. Since 
the isothermal compressibility kt is indeed a thermodynamic manifestation of the fluc- 
tuations in the particle density, isotherms p, = p{T,p) contain fundamental information 
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on the nature of the phase transition involved [TB]. 

From a more general point of view the equation of state /i = yu(p, T) of an interact- 
ing quantum gas contained in a rigid vessel of volume V has a fundamental relevance on 
itself. Such a quantity enables the prediction of the density profile n{f) of the interacting 
gas confined by an inhomogenous trap potential V^xti^ which is the main measurable 
property in the current experiments with ultracold alkali gases[71 |Hl [9], |T7j. The connec- 
tion between the thermodynamics of the gas in a box of volume V and the thermody- 
namics of the gas confined within an external potential Vext (^) can be achieved through 
the Local Density Approximation (LDA). In such procedure one obtains p = p{fi,T) 
by inversion of the equation of state, and then replacing fi hj fi — Ve^t(r). This gives 
the density profile n{r) = p{fi — Vextif), T). This procedure has been shown to be exact 
in the appropriate thermodynamic limit of a gas confined by a trap Ka;t(?^[IHl [iHl 120] . 
As we shall carefully discuss, if the trapped gas is in the BEC phase the density profile 
bears the information of the isotherm p(//, T) above, below and at the transition. Thus, 
knowledge of yU = fi{p,T) also yields the equation of state for the trapped gases. Inac- 
curate calculations of the equation of state in the homogenous case will be inherited to 
the inhomogeneous ones. One of the purposes of this article is to suggest high-resolution 
measurements of the density profiles in the current experiments. These would not only 
settle the issue of the validity or not of the mean-field calculations but would also pave 
the way to improve theoretical descriptions and thus better understanding of the BEC 
phase. 

The paper is thus organized as follows. In Section II a general discussion of the ther- 
modynamics of phase transitions and the weakly interacting gas is presented. Section 

III is devoted to a brief presentation of the approximations considered in this work and 
the results for the equation of state /i = /i(T, p) corresponding to each of them. Section 

IV discusses the implications of these results within the general theory of the weakly 
interacting gas and in the interpretation of current experiments in ultracold alkali gases. 

2. General thermodynamic considerations and the weakly interacting Bose 
gas 

For our analysis, we consider the equation of state p = fi(T,p), that is, the dependence 
of the chemical potential p, on temperature T and particle density p = N/V, with N 
the number of particles and V the volume of the sample. The extensivity property 
of Helmholtz free energy F = F{N,V,T) allows us to either write F = Vf{p,T) or 
F = Nf{v,T) with V = V/N = p~^ such that f = pf, and both forms carry the same 
physical information. One can either calculate the chemical potential as 




(1) 



Order of BEC in mean-field approximations 



4 



or the pressure 

(2) 

Therefore, /i = /i(T, p) and p = p{v, T) also carry the same information and the phase 
diagram may be inferred from any of these forms. The laws of thermodynamics require 
that both of those expression are single valued, that is, for any given values of p and T, 
or V and T, there must only exist a single value of or p respectively. In addition, these 
formulae carry out crucial information on the order of the phase transition involved. The 
stability of a particular thermodynamic state is ensured if the isothermal compressibility 
kt is positive, where kt may be either calculated as 



or 

-1 




(3) 
(4) 

These equations imply that {dp/dp)T > and {dp/dv)T < for any thermodynamic 
equilibrium state. An exception is a critical point where the latter derivatives become 
zero. 



The emergence of a negative region for the isothermal compressibility signals the 
onset of a first-order phase transition where two phases with different densities, entropies 
and energies coexist at the same temperature, pressure and chemical potential [IE]. Let 
us consider here that such states have densities pi and p2- The equal- areas Maxwell 
construction enables us to calculate the values of the thermodynamic properties of 
interest for both coexisting states. In order to see this, it is enough to consider the 
Gibbs-Duhem relationship: 

- SdT + Vdp - Ndp = 0, (5) 

where S is the entropy. By assuming the system is analyzed along an isotherm, we can 
set dT = and look for the coexisting values of the density or the volume. Then, we 
can either impose equal chemical potential p = Pcoex at both phases yielding 



J vdp = j dp = 



V2 



Pcoex{v2 -Vi) - / pdv = 0, (6) 



where Pcoex, vi and V2 are the volume values at coexistency, with vi = p^^ and V2 = P2^. 
The value of the chemical potential is Pcoex = /^(^)Pi) = l^iT,p2). Equation ([6]) is 
the usual Maxwell equal-area construction in the p — v diagram. Alternatively, if one 
imposes equal pressures at coexistence, namely p = Pcoex, one finds from Eq.([5]) 

J pdp = J dp = 

= l^coex{p2 - Pi) - / pdp = 0, (7) 
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with Hcoexi Pi and p2 the corresponding values at coexistency. Equation ([7]) is also a 
Maxwell equal-areas construction but in the fi — p diagram. 

The onset of a second-order phase transition is featured by the divergence of 
the isothermal compressibility kt and the heat capacity Cy characterized by universal 
exponents|16[[2T|[22]. In this work we focus our attention to the isothermal compressibil- 
ity whose divergence can be expressed by the vanishing derivative {dfi/dp)T = or 
{dp/dv)T = 0, at the transition. This condition demands that the isotherm p = p(T,p) 
must become "flat" at the critical density p = pc- We shall explicitly show that none 
of the analyzed approximations for the weakly interacting Bose gas show this strong 
requirement. It is worth recalling that an ideal Bose gas confined in a box of volume V 
does show a diverging compressibility at the critical temperature T^pT| [23]. 

The theory of the weakly interacting Bose gas assumes that the atoms have no 
structure since they typically are in the same hyperfine state |9J. If a pairwise interatomic 
potential U = ?7(|rj — r^|) is assumed the Hamiltonian in second quantization can be 
written as 

H = Y1 4^l^k + ^ E ^ iq)al+gal_^akak', (8) 

where k is the three dimensional wavevector quantized in a box of volume V, a|, and 
are creation and annihilation operators of particles with momentum hk, e° = U^k"^ /2m 
is the one-particle kinetic energy, and U{q) is the Fourier transform of the interparticle 
potential U{r). The main assumption for the description of an interacting gas at low 
temperatures is the contact potential approximation, namely, U{q) ~ Uq for all g, with 
Uq = Anfi^a/m representing the strength of the interaction and a the s-wave scattering 
length which in this work we shall assume as positive. In the t-matrix approximation 
also considered here, corrections to this simple potential are also included [2]. An ad- 
ditional ansatz completes the framework. A Bose gas is considered weakly interacting 
if pa^ <^ 1. All the approximations considered in this work fgall satisfy these general 
requirements. In the following we succinctly present the differences between each of 
them that yield to different equations of state /i = /i(p, T). We refer to the reader to 
the original sources cited along this work for further details. 

For our calculations we use units h = m = a = 1. We shall analyze three isoterms 
in HF, P and YY approximations, ksT = 0.01, 0.1 and 1.1 in dimensionless units, while 
only ksT = 0.1 in TM approach. Before proceeding to the calculations it is worth to 
state some words on the physical regime associated to the isotherms chosen. Typical 
experiments on ^^Na ultracold gases [17] show that BEC transition temperature is nearly 
Tc ~ 100 nK with a scattering length a ^ 55ao, where ao is Bohr radius. This yields 
in dimensionless units kBT ^ 10~^, thus indicating that the chosen temperatures are 
a bit our of range, but this is not necessarily so. As we shall discuss below, while in 
HF and P the behavior is qualitatively the same for all temperatures, YY does show 
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three different qualitative beliaviors corresponding approximately to those temperatures. 
Moreover, recent advances in the experimental achievement of ultracold atomic gases 
have shown that the scattering length a can be tuned by external magnetic fields to 
larger values near a Feshbach resonance [211 ISSl [261 127]. Hence, if an increase of tenfold 
in the scattering length is achieved, then a temperature of T f» 100 nK would give fc^T ^ 
0.01 in dimensionless units which corresponds perfectly to the first isotherm considered 
in this work. The tunability of the scattering length a together with with the fact that 
each time is possible to achieve condensates with a larger number of particles suggest 
that in the near future BEC in ultracold gases will be achieved for higher temperatures. 
Its is important to underline that the criterion for a weakly interacting gas is fulfilled 
here. The gas parameter associated to the isotherms are pa^ = 0.00017, 0.0052 and 
0.19, respectively. Admittedly, the last value being at the border of what one should 
consider the g weakly interacting. 

3. BEC transition within four mean-field theories 

In this section we analyze four different approximations that describe the weakly inter- 
acting Bose gas. These schemes have provided accurate descriptions at and near zero 
temperature. Given this success their use have been extended to finite temperatures near 
criticality. As suggested before, BEC can be attained either at fixed density lowering 
the temperature or at fixed temperature increasing the density. In this work, we shall 
make use of the latter approach and thus BEC transition will be accomplished in /i — p 
space at fixed temperature T by varying the total density. In this case, the transition oc- 
curs at a critical density pc, the gas being normal for p < Pc while BEC sets in for p > pc. 



3.1. Hartree-Fock 

The Hartree-Fock approximation is a self-consistent approach in which the state of the 
A^— particle system is expressed in terms of effective one-particle states, yielding a gas 
of non-interacting excitations whose energy spectrum depends self-consistently on its 
density and the actual interparticle interaction [101 111]- Then, one can construct the 
grand potential as = (H) — TS + yu(A^), in terms of variational occupation numbers 
/q, with 

(^) - ^ + E {4 + 2pUo) fk-^T. hu (9) 

and 

{N)^No + J2h (10) 

with A"o the number of particles in the condensate with k = 0. Further, one assumes 
that the entropy is given as that of an ideal Bose gas, but in terms of the effective 
one-particle occupation states fk, this is S' = kJ2k [(1 + fk) ln(l + fk) — /fcln/fc]. At 
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zero temperature, HF is considered an adequate approximation being equivalent to the 
Gross-Pitaevskii approach. However, at low but finite temperatures it predicts an exci- 
tation spectrum with a gap = e° + 2p[/o as seen from Eq.(l9]) contrary to the expected 
gapless spectrum [28j. that P and YY approximations do take into account. However, 
in the normal gas region, that is, for densities below and at the transition, HF is an 
acceptable approximation. 

By looking for the values of fk that minimize Vt at constant temperature, one is 
lead to the thermodynamic equations in HF approximation. One finds that for densities 
p < Pc the gas is in the normal phase with|10l [TT] 

p=^g3/2[Hp-'^pUo)]- (11) 

Anp 

On the other hand, for p > Pc BEC sets in, and it is found that, 

P = Po + T393/2 {-PpqUq) with /i = (2p - po) Uo, (12) 
where p and po are the total and condensate densities. In the above equations 



At = h/y/2TTmkBT, and we used the Bose integral 

9n{a) = / -dx, 13 

1 [n) Jo ° — 1 

with P = l/keT. 

The onset of BEC is assumed to occur at po = in Eq.(fT2l). namely at pc = 2pcf/o. 
This defines a relation between T and pc, which is the same as in the ideal gas, 

Pc = ^93/2{0). (14) 
Arp 

Eqs. ffTTl) and f|T2|) yield p = fi{p,T) as a continuous function of p at BEC. As we now 
show, this continuity does not guarantee the existence of a second-order phase transition. 

Fig. [Ushows solutions to Eqs. (ITT]) and (fT2!) for the three isotherms, ksT = 0.01, 0.1 
and 1.1. The behavior of p vs p for all temperatures is essentially the same. In the normal 
gas region, p < pc, the derivative {dfi/dp)T 2f/o as p — )■ pc. This straightforwardly 
predicts that the isothermal compressibility does not diverge at BEC. The behavior 
of p vs p corresponding to the BEC region, i.e. the solution to Eq.( [T2|) . shows an 
anomalous thermodynamic behavior: this solution, in principle only valid for densities 
p > Pc, intrudes into the normal region yielding a multiple valued chemical potential 
as a function of density. Thermodynamics forbids this multiple-valuedness. However, 
as done originally in Refs.|I^ and [I5], such a behavior can be avoided by invoking 
a first-order phase transition, joining two phases with different densities by means of 
an equal-areas Maxwell construction, as shown in Fig. [T]for the isotherm = 1.1. 
As discussed above, equal-areas Maxwell construction ensures that the two phases with 
different densities have the same temperature, pressure and chemical potential, namely. 
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Figure 1. Isotherms /i vs p at fc^T = 0.01, 0.1 and 1.1, from numerical solutions of HF 
equations (|lll) and (|12l) . /ic and pc are the values at BEC. The gaseousness parameter 
at criticality are pcO? — 0.00017, 0.0052 and 0.19, respectively. The isotherm fc^T = 1.1 
shows Maxwell equal-areas construction for a first-order transition. 

that the phases coexist. This is the signature of a first order phase transition and one 
finds that the conditions for a second-order phase transition, as discussed before are 
never met. As described previously the fact of having the coexistence of two states with 
different densities imphes two different entropies, namely a latent heat, and two values 
of the energy. In the same fashion, a multiple valued chemical potential gives rise to 
also an unacceptable multiple- valued pressure of the system. 

3.2. Popov 

The Popov approximation is considered correct at zero temperatures being essentially 
the same as Bogoliubov approximation [28] . This approximation, as shown below, yields 
a gapless excitation spectrum [21 [31 ITTj . linear in the excitation momentum as — )■ 0, and 
that gives rise to the phenomenon of superfluidity, just as in ^He[6]. Bragg spectroscopy 
in an ultracold ^'^Rb gas[29j has shown that the excitation spectrum is clearly given by 
the Bogoliubov expression, thus providing validity to this approximation at very low 
temperatures and, incidentally, confirming the superfiuid nature of these gases at the 
BEC state. 
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The number- conserving Hamiltonian of the system at this level of approximation 
is given by[TT] 

E (4aLfc + afcfl-fc) . (15) 

It looks very similar to the Hartree-Fock version except for the last term which represents 
the annihilation of two particles into the condensate and the creation of two particles 
from the condensate. For a finite temperature calculation of the thermodynamics once 
again one deals with the problem in the grand canonical ensemble. After performing 
the usual Bogoliubov transformation pTj. the Popov procedure follows a similar line as 
HF and the system emerges as a kind of ideal gas with elementary excitations whose 
excitation spectrum is, 



ek = ^ ^ + 2pUo-fi -pWI (16) 

At a fixed temperature T and for densities below the critical one pc, the equation in 
P approximation for the density p as a function of temperature T and chemical potential 
/i, turns out to be the same as in HF, Eq.( ITT|) . but for densities p > Pc it is found that 
the density of the fluid is 

2 {x + (3poUo) x^/^dx 

^ ^° ny^X^Jo y/x^TWMJox ^^^^^+WJ^ _ I ^ ' 

while the chemical potential is again given by /x = (2p — po)Uo as in HF. The transition 
occurs when the condensate density vanishes, po = 0, and the critical density pc is 
given by the condition f lT^ . The solution n{p,T) for p > pc becomes continuous at 
Pc with the solution to Eq.lITT]). Fig. [2] shows three isotherms p vs p obtained fo the 
Popov approximation using Eqs. (fTT!) and (ITTj) for the same values used before. It is 
found that although the multiple-valued region is smaller in Popov description than 
in Hartree-Fock the thermodjTiamic behavior is qualitatively similar. An equal-areas 
Maxwell construction would yield to a flrst-order phase transition as well. 

3. 3. Yukalov- Yukalova 

This recent scheme is based on the inclusion of anomalous averages Ui = {^i{r)-ipi{r)) 
into the Hamiltonian. According to its authors the anomalous averages should not be 
neglected since their contribution is of the same order as that of the condensate density 
[T2] . The operator ^/'i(r) is given in terms of the creation operators as 

i'i{^ = Y.MrW (18) 

fcT^O 

It is discussed in Ref. [12] that these anomalous are not properly considered in Popov 
approximation and therefore their inclusion provides with a better description for the 
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Figure 2. Isotherms /i vs p at ksT = 0.01, 0.1 and 1.1, from nmnerical solution of 
Popov equations ([TT|) and p7)) . and pc are the values at BEC. The gaseousness 
parameter at criticality are pcO? — 0.00017, 0.0052 and 0.19, respectively. 

BEC transition. In order to obtain the thermodynamics of the system YY approxima- 
tion considers exphcitly the symmetry breaking of the state by the appearance of the 
condensate and implements a so-called representative ensemble calculation where one 
Legendre multiplier is associated to the condensate particles po and a different one to 
the non-condensate fraction p — po- This takes into account correlations of uncondensed 
particles and allows for a clear distinction of the superfiuid and condensate fractions 
at finite temperatures. At zero temperature, this approximation also agrees with that 
of Bogoliubov|28] as expected. The details of the calculation can be consulted in the 
original works of Yukalov and Yukalova|12]. 



Once more, nevertheless, for a given temperature T and in the normal region p < pc, 
the relevant equations are those of HF, Eq.f jTTj) . For densities above criticality, p > pc, 
a new set of equations are given for the the value of the density p and the anomalous 
average ai namely. 



P 




Q poo . , 



1/2 



dx 



(19) 
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mc\'' 



PqUq 



x^ 



1/2 



1/2 



X 



x^ 



coth 



mc^x' 



2kBT 



1 



dx 



(20) 



where the speed of the elementary excitations c is given in terms of the condensate 
density po and ai by mc^ = (po + cti) Uq. In this approximation, the chemical potential 
is given by [12] 

p/i = Po (2p - Po + cTi) f/o + (p - Po) (2p - Po - cTi) f^o- (21) 

In this case BEC sets in when both ai and po vanish. This yields pc = 2pcf/o at the 
transition with p^ determined by the condition f|T^ . As in HF and P, the solution p(p, T) 
to ffT^ - fl^Tl) is continuous with the solution of flTT]) at BEC. This approximation is much 
richer than the two previous ones. The behavior of the isotherms p vs p show three 
different regimes depending on the temperature, see Fig. [31 At very low temperatures, 
i.e. ksT = 0.01, there are no multiple valued solutions, neither an unstable region. 
Thus, there is no need for appealing to a first-order phase transition. However, despite 
the fact that the transition appears continuous, it cannot be considered as a bona-fide 
second order one since neither side shows a divergent compressibility at BEC; rather, 
the compressibility is discontinuous. Nevertheless, this case of YY is the closest to a 
second-order phase transition. At intermediate temperatures, i.e. ksT = 0.1, although 
no multiple valued solutions exist, an unstable region appears where {dfi/dp)T < 0, 
and a first-order phase transition must therefore be adscribed. At large temperatures, 
/c^T = 1.1, the multiple- value issues of HF and P pervade YY as well. 



3.4- Many-Body t-matrix 

This approach was originally carried out in Ref. [13] but here we follow the analysis 
given in Ref. [2]. This approximation builds on the Popov theory and takes into ac- 
count higher corrections to the scattering t-matrix including many-body effects, thus 
going beyond the contact interaction potential U{k) ^ Uo- 



This case does not have HF as the solution in the normal phase. Instead, for p < Pc 
and fixed T, the density equation for the normal gas is, 

P=^^73/2(/3A), (22) 



A 



T 



while for p > Pc 



P = Po + ^ 



2 



X 



dx. 



In the above equations, A is determined by an additional quantity. 



E.. 1 



(23) 
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Yukalov-Yukalova Approximation 




Figure 3. Isotherms ^ vs p at ksT = 0.01, 0.1 and 1.1, from numerical solution of 
YY equations PT|) and p^ - (pij) . /ic and Pc are the values at BEC. The gaseousness 
parameter at criticality are pcU^ — 0.00017, 0.0052 and 0.19, respectively. 



where, for p < pc, = x — /3A, and for p > p^, E^ = a/x^ + 2/3Aa;. These 
equations must be solved self-consistently with the corresponding equation for the 
chemical potential, 

2pUo 
(1 + alio) 

BEC occurs when A — and a — )■ oo from both sides. This yields /i = at the transi- 
tion, as in the ideal gas, and the transition density is again given by condition (IT4|) . 



Figure m shows the isotherm ksT = 0.1, with pca^ = 0.0052. Given the fact that 
this theory goes beyond mean-field by considering the renormalization of the coupling 
parameter, it is somewhat surprising that its predictions, in a way, fare less satisfactory 
than HP and P. It is seen that in the normal region p < pc the solution is unstable near 
BEC, i.e. {dfi/dp)T < 0. In the BEC side, p > Pc the situation is more worrisome, since 
it appears that {dfj,/dp)T becomes extremely large, probably diverging at the transition. 
This indicates that the compressibility becomes nearly zero, that is, the gas becomes 
incompressible at BEC. This is completely opposite to a critical behavior where large 
density fluctuations are due to the large compressibility of the gas [22]. 
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Figure 4. Isotherm fivs p at ksT = 0.1, from numerical solution oit matrix equations 
At BEC this approximation predicts a full incompressibility of the gas. This 
property is completely unexpected and it is contrary to the standard results in weakly 
interacting systems. 

4. Final remarks and Perspectives 

In this work we considered four of tlie most used mean field tlieories to describe the 
onset of the BEC phase in weakly interacting gases. We calculated the equation of 
state /i = /i(p, T) corresponding to each of these theories at three different tempera- 
ture regimes (with exception of the t-matrix theory which has calculated only at one 
temperature) satisfying the requirements of dilute gas in the contact potential approx- 
imation. The central conclusion is that none of the descriptions considered display a 
true second-order phase transition at BEC as expected by general symmetry break- 
ing considerations. In addition and excepting the Yukalov-Yukalova approximation at 
ksT = 0.01 and 0.1 all the approximations considered exhibit an unphysical behavior 
in the vicinity of BEC. This anomalous property restricts severely its applicability at 
finite temperatures near the transition. 

The interest in the equation of state equation of state /i = /i(p, T) goes beyond the 
order of the phase transition involved. As discussed in the Introduction it provides a 
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fundamental link to the determination of the local density profile n{r) in gases confined 
by inhomogeneous magnetic or optical traps which is one of the main measurable quan- 
tities in the current experiments with ultracold alkali gases[7l El |30]. It is worthwhile to 
point out that a single profile at a temperature below BEC in the trapped gas bears in- 
formation of the homogeneous one for densities p below, at and above the critical one pc- 

Using the Local Density Approximation (LDA) we have calculated density profiles 
for an isotropic 3D harmonic trap, Vextif) = \muj'^r'^ using the isotherms p = p{p,T) 
obtained from the analyzed theories. In this case the thermodynamic limit corresponds 
to — 7- oo, a; — 7- but Nu^ = constant. It is found that the predicted density pro- 
files for HP and P are multiple valued for any temperature, while for YY this occurs 
only for relatively high temperatures. As an illustration of these multiple valuedness 
we show in Fig. [5] the YY profile obtained for ksT = 1.1 by applying LDA to the cor- 
responding case in Fig. [31 We recall that LDA, within the HF and P approximations, 
has been widely used to calculate thermodynamic properties of trapped gases, see e.g. 
Refs.[9l Uni inn [321 [33], and the issue of multiple- valuedness has not been brought up. 
This is perhaps due to the fact that, when performing numerical calculations, the mul- 
tiple valued region being very small can be "short-circuited", inadvertently or not. In 
that case, the resulting profile is thus continuous and, apparently, the ensuing thermo- 
dynamic properties of the trapped gas are quite insensitive to this correction. We expect 
that the results provided in this work will bring again the attention to this matter. An 
accurate determination of density profiles would help to settle this issue and certainly 
it would provide solid guidelines to a better formulation of the theory of the weakly 
interacting gas. 

Although the mean-field theories here described do not predict a second-order phase 
transition, their behavior near T = remains correct and valuable. With the exception 
of HF, the others correctly incorporate the gapless elementary excitation spectrum, 
originally predicted by Bogoliubov[28]. The results presented here suggest strongly that 
at temperatures near the transition additional properties of the interatomic potential 
are needed to fully capture the thermodynamics of the system. Nevertheless it is 
important to stress that exact calculations of the thermodynamics of quantum systems 
based on complete solutions of the Hamiltonian ([8]) are beyond present capabilities. We 
believe that in the construction of a new approximation or an improvement of those, the 
calculation of the equation of state p = p{p, T) should be a useful tool to determine the 
predicted order of the transition, and therefore, of the correct thermodynamic behavior 
of the ultracold gas in the BEC state. 
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Figure 5. Density profile p{r) of a harmonically trapped Bosc gas in the YY 
approximation at temperature UbT = 1.1 obtained through the implementation of 
Local Density Approximation. 
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